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Abstract 

The discovery of Massive Astrophysical Compact Halo Objects (MA- 
CHOs) in microlensing experiments makes it compelling to understand 
their physical nature, as well as their formation mechanism. Within the 
present uncertainties, brown dwarfs are a viable candidate for MACHOs, 
and the present paper deals with this option. According to a recently 
proposed scenario, brown dwarfs are clumped along with cold molecular 
clouds into dark clusters - in several respects similar to globular clusters 
- which form in the outer part of the galactic halo. Here, we analyze 
the dynamics of these dark clusters and we address the possibility that a 
sizable fraction of MACHOs can be binary brown dwarfs. Moreover, we 
point out that Ly-a absorption systems naturally fit within the present 
picture. 

Subject headings: Dark matter - Galaxy: halo - Galaxy: kinematics and 
dynamics - Gravitational microlensing 
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1 Introduction 



Observations of microlensing events (Alcock et al. [ 1993] , Aubourg et al. |1993|) 



towards the Large Magellanic Cloud (LMC) strongly suggest that a substantial 
fraction of the galactic halo should be in the form of dark compact objects, 
named MACHOs (massive astrophysical compact halo objects) (Dc Riijula, Jet- 
zer and Masso |l992fl ). 



Actually, the MACHO collaboration has recently announced the discovery 
of several new events during their second year of observations (Alcock et al. 
[1997Q, and 8 microlensing events have been detected so far []. 



Although the presently-available limited statistics prevents clear-cut conclu- 
sions to be drawn from experimental data (Gates, Gyuk and Turner [ 199(| ). 



the evidence for such a discovery is firm and its implications are striking. In 
fact - under the assumption that MACHOs are indeed located in the galac- 
tic halo - the inferred halo mass in MACHOs within 50 kpc turns out to be 
2.01q y x 10 11 M© (Alcock et al. [1997]), which is several times larger than the 
mass of all known stellar components of the Galaxy and represents a relevant 
portion of the galactic dark matter. Remarkably enough, this result is almost 
independent of the assumed galactic model. Unfortunately, this circumstance 
contrasts with the strong model-dependence of the average MACHO mass. It 
has become customary to take the standard spherical halo model as a baseline for 
comparison. Regretfully - because of the low statistics - different data-analysis 
procedures lead to results which are only marginally consistent. Specifically, 
within the standard halo model, the average MACHO mass reported by the 

1 It should be mentioned that the MACHO team has found at least seven more events 
(which are reported on the Alert list), but a full analysis of them has not yet been published. 
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MACHO team is 0. 461^2 M Q (Alcock et al. [1997]), whereas the mass moment 
method (De Rujula, Jetzer and Masso [1991]) yields 0.27 M Q (Jetzer [1996]). 

What can be reliably concluded from the existing data-set is that MACHOs 
should lie in the mass range 0.05 — 1.0 Mq (see also Table 9 in Alcock et al. 
[1997]), but stronger claims are unwarranted because of the high-sensitivity of 
the average MACHO mass on the uncertain properties of the considered galactic 
model (Evans [1996], De Paolis, Ingrosso and Jetzer [1996]). 

Mass values > 0.1 Mq suggest that MACHOs should be either M-dwarfs or 
else white dwarfs. Observe that these mass values naturally arise within the 
standard halo model. 

As a matter of fact, deeper considerations show that the M-dwarf option can 
look problematic. The null results of several searches for low-mass stars both 



in the disk and in the halo of our galaxy (Hu et al. [1994]) entail that the halo 
cannot be mostly in the form of hydrogen-burning main-sequence M-dwarfs. 
Optical imaging of high-latitude fields taken with the Wide Field Camera of 
the Hubble Space Telescope indicates that less than ~ 6 % of the halo mass 



can be in this form (Bahcall et al. [ 1994 1 ) . A more detailed analysis which 
accounts for the fact that halo stars likely have a lower metallicity (with respect 
to the solar one) leads to an even more stringent upper limit of less than ~ 1 % 



(Graff and Freese [ 1996 ). We emphasize that these results are derived under 
the assumption of a smooth spatial distribution of M-dwarfs, and become less 



severe in the case of a clumpy distribution (Kerins [ 1996 1 ) . In the latter case, 



as pointed out by Kerins [ 1997 1 , the dynamical limits and HST observations 
require that the overwhelming fraction of M-dwarfs, at least 95%, must still 
reside in clusters at present. His analysis shows that there exists a wide range 
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of cluster masses and radii which are consistent with these requirements. 

As we said, an alternative explanation for MACHOs can be provided - 
within the standard spherical halo model - by white dwarfs, and a scenario 
with white dwarfs as a major constituent of the galactic halo dark matter has 
been explored (Tamanaha et al. [ 1990 1 , Fields et al. [ 1996| , Adams and Laugh- 



lin 1996 , Chabrier, Segretain and Mera [ 1996 1 ) . However, even this proposal 
meets difficulties. Besides requiring a rather ad hoc initial mass function (IMF) 
of the progenitor stars sharply peaked somewhere in the range 1 — 8 Mq and 
a halo age larger than ~ 16 Gyr, strong constraints on the number density of 
halo white dwarfs arise from present-day metal abundances in the interstellar 



medium (Ryu, Olive and Silk [199C], Gibson and Mould [1996]) and from deep 
galaxy counts (Chariot and Silk | 1995fl ). In any case, future HST deep field 
exposures will either find the white dwarfs or put constraints on their fraction 
in the halo (Kawaler J1996 ). 

Mass values <0.1 Mq make brown dwarfs ^] an attractive candidate for 
MACHOs 0. In fact, these mass values are supported by several nonstandard 
halo models. An example is provided by the maximal disk model (van Albada & 
Sancisi [1986]; Persic & Salucci [1990]; Sackett [1996]), in which more matter is 
comprised within the disk whereas the halo is less massive as compared with the 
standard halo model. The latter fact implies a falling rotation curve, and so a 
smaller transverse velocity of MACHOs. Hence, the microlensing time scale gets 
longer for a given MACHO mass, which means a smaller implied MACHO mass 

2 Brown dwarfs have been disc overe d quite recently in t he sol ar neighbourhood and in the 
Pleiades cluster (Rebolo et al. |l995| , Nakajima et al. [|199§). The idea that MACHOs 
are brown dwarfs has been contemplated by several authors (see Carr [ 1994 and references 
therein). 

3 We notice that the limit for hydrogen burning — usually quoted as 0.08 Mq - gets in- 
creased up to 0.11 Mq for low-metallicity objects, such as a halo population (D'Antona [1987], 
Burrows, Hubbard & Lunine [1989]). 
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for a given observed timescale. We stress that a reduced transverse velocity of 
MACHOs arises also in other nonstandard halo models. For instance, a radially 
anisotropic velocity distribution or a halo rotation would do the job (Alcock et 
al. [1997]; Evans [1996]; De Paolis, Ingrosso & Jetzer [1996]). We also notice 



that the EROS collaboration (Renault et al. |1997fl ) has shown that MACHOs 
in the mass range 10 -7 — 2 x 10~ 2 M Q do not contribute significantly (less 
than 20%) to the halo dark matter (this result is consistent with the MACHO 
experiment for objects of mass O.IMq to IMq). 

Although present uncertainties do not permit to make any sharp statement 
about the nature of MACHOs, brown dwarfs still look as a viable possibility to 
date, and we shall stick to it throughout. 

Even if MACHOs are indeed brown dwarfs, the problem nevertheless remains 
to explain their formation, as well as the nature of the remaining dark matter 
in galactic halos. 

We have previously proposed a scenario in which dark clusters of brown 
dwarfs and cold molecular clouds - mainly of H2 - naturally form in the halo 
at galactocentric distances larger than 10 — 20 kpc (De Paolis et al. |l995a - 



[ 1995d ] ) . Similar ideas have also been put forward by Gerhard and Silk [ 1996 ] . 
A slightly different picture based on the presence of a strong cooling-flow phase 
during the formation of our galaxy has been considered by Fabian and Nulsen 



[1994, 1997 1 and leads to a halo made of low-mass objects. In addition, Pfen- 



niger, Combes and Martinet [ 1994 1 suggested that H2 clouds can constitute the 
dark matter in the disk of our galaxy. 

The model in question encompasses the one first proposed by Fall and Rees 



[1985 to explain the formation of globular clusters, and no substantial additional 
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hypothesis is required. Various resulting observational implications have also 
been addressed. In particular 

i) the 7-ray flux arising from halo molecular clouds through the interaction with 



high-energy cosmic-ray protons has been estimated (De Paolis et al. [ 1995a , 



1995J); 

ii) an anisotropy in the Cosmic Background Radiation (CBR) is predicted to 



show up when looking at the halo of M31 galaxy (De Paolis et al. [ 1995c ]); 
Hi) the infrared emission from MACHOs located in the halo of M31 galaxy 
turns out to be observable with the detector on ISO satellite or with the next 



generation of satellite-borne detectors (De Paolis et al. [ 1995c |). 

We would like to stress that a large amount of MACHOs - up to 50% in 
mass - can well consist of binary brown dwarfs, formed either by the same 
fragmentation process that produces individual brown dwarfs or later when the 
dark clusters start to undergo core collapse. 

The aim of the present paper is to discuss in a systematic fashion further 
aspects of the above scenario. Basically, we try to figure out the dynamics of 
dark clusters. More specifically, we investigate the constraints which ensure 
their survival against various kinds of gravitational perturbations. In addition, 
we demonstrate that - because of dynamical friction on molecular clouds in 
the dark cluster cores (to be referred to as frictional hardening) - the present 
orbital radius of not too hard binary brown dwarfs turns out to be typically 
of the order of the Einstein radius for microlensing towards the LMC. As a 
consequence, taking also into account the adopted selection procedure in the 
data analysis, we understand why they have not been resolved so far - still, 
we argue that they can be resolved in future microlensing experiments with a 
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more accurate photometric observation. Finally, we show that hy-a absorption 
systems naturally fit within our model. 

The plan of the paper is as follows. In Section 2 we recall the main points 
of the considered picture for the formation of dark clusters. Various dynamical 
constraints are thoroughly analyzed in Section 3, paying particular attention to 
the phenomenon of core collapse. In Section 4 we study the process whereby 
brown dwarfs form close binary systems (as a consequence of core collapse) and 
we investigate the mechanism of frictional hardening in great detail. In Section 
5 we turn our attention to the thermal balance in the halo molecular clouds. 
Section 6 contains a short discussion of the relevance of Ly-a absorption systems 
for the present scenario. Our conclusions are offered in Section 7. 

2 Scenario for dark cluster formation 



As shown elsewhere (De Paolis et al. [1995a, 1995t]), the model in question 



encompasses the one first considered by Fall and Rees [1985] for the formation 



of globular clusters and relies upon the result of Palla, Salpeter and Stahler 



[ 1983 (hereafter PSS) that the lower bound on the Jeans mass in a collapsing 
metal poor cloud can be as low as 1O -2 M0, provided certain environmental 
conditions are met. 



Let us begin by summarizing the ideas of Fall and Rees [ 1985 1 from a point 
of view which is most convenient for our considerations. 

After the initial collapse, the proto galaxy (PG) is expected to reach a quasi- 
hydrostatic equilibrium state with virial temperature ~ 10 6 K. Fall and Rees 



4 A so mewh at different extension of the Fall and Rees [ 1985 scheme has been proposed by 
Ashman JL99q| . Amon g h s motivations was the large spread in the age of globular clusters 



found by 



Larson [119901 . However, it seems nowadays that the s 



have a small scatter in ages (Unavane, Wyse and Gilmore [ 1996 1 ) 



ars in the halo of our galaxy 



G 



[1985 1 have shown that in such a situation a thermal instability develops: density 
enhancements rapidly grow as the gas cools to lower temperatures. In fact, 
irregularities in the inflow during the gas collapse and also fluctuations in the 
distribution of nonbaryonic dark matter (if present on the galactic scale) would 
introduce perturbations with a wide range in size and amplitude. As a result, 
randomly distributed overdense regions will form inside the PG. For reasons 
that will become clear later, these overdense regions will be referred to as proto 
globular cluster (PGC) clouds. 

Under the assumption that the plasma in the PG is in collisional ionization 
equilibrium, it turns out that the cooling rate (as a function of density p and 
temperature T) has the form 

A(p,T) = p 2 L(T) (1) 



(the expression of L(T) can be found e.g. in Einaudi and Ferrara [1991]). The 
cooling time is 

_ 3pk B T 

cooi ~ 2p(A - T) ' 1 ' 

where the heating rate T due to external heating sources has been taken into 
account (here p ~ 1.22 m p is the mean molecular mass of the primordial gas). 
Since at the high temperatures under consideration the heating rate can safely 
be neglected, it follows that t coo i ~ p . On the other hand, the free-fall time 
reads 

t f/ ~ {Gp)-^ , (3) 

and so we see that t coo i decreases faster than tff as p increases. As the above 
quasi-hydrostatic equilibrium state of the PG is characterized by the condition 
tcooi ~ ttf, it is clear that inside the PGC clouds we have t coo i < tff. That is 
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to say, the PGC clouds cool more rapidly than the rest of the PG. This process 
continues until hydrogen recombination occurs, because as soon as this happens 
- at a temperature ~ 10 4 K - the cooling rate decreases precipitously, under the 



ionization equilibrium assumption (Dalgarno and McCray | 1972 |). Therefore, 
the regime t coo i > t// should now be established even in the PGC clouds and so 
the PG can be regarded at this stage as a two-phase medium, with cold PGC 
clouds in pressure equilibrium with the external (inter PGC clouds) diffuse hot 
gas. 



However, Kang et al. 1990 realized that the fast radiative cooling of the 
PGC clouds (from 10 6 K to 10 4 K) implies that the plasma inside these clouds 
cools more rapidly than it recombines, so that the above ionization equilibrium 
assumption is violated. Actually, the out-of-equilibrium recombination results 
in an enhanced ionization fraction. This fact does not affect the previous con- 
clusions for temperatures > 10 4 K, but entails drastic changes at lower temper- 
atures. Indeed, the existence of a sizable amount of protons and electrons at 
temperatures < 10 4 K gives rise to H 2 formation via the following reactions 

H+p^H++ ll H + e^H-+ 1 (4) 

and 

H+ + H^H 2 +p , H + H- -» H 2 + e . (5) 

As a consequence, the PGC clouds undergo a further cooling below 10 4 K. 
Specifically, there is a direct radiative cooling via reactions (|J) and a radiative 
cooling via excitation of roto-vibrational transitions of H 2 (observe that H 2 is 
further produced by reactions (||) and @ below). We stress that the latter 

5 We wish to emphasize that the familiar reaction H + H — » H2 + 7 (which requires grains 
in order to be efficient) is presently irrelevant, since no dust exists in the metal poor PGC 
clouds. 
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process is very effective (much more than the former) at temperatures < 10 4 K 
and plays a crucial role in our considerations. Since now t coo i <C tft in the PGC 
clouds, the collapse goes on and the PGC cloud density rises steadily. When 
the number density in the PGC clouds exceeds 10 s cm -3 , the H 2 production 
increases dramatically thanks to the three-body reactions 

H + H + H^H 2 + H, H + H + H 2 ^H 2 + H 2 , (6) 

as pointed out by PSS. In effect, these reactions are so efficient that virtually 
all the atomic hydrogen gets rapidly converted [] to H 2 . Correspondingly, the 
cooling of the PGC clouds is strongly enhanced and their evolution can proceed 
as in the scenario proposed by PSS. 

Still, it goes without saying that H 2 can be dissociated by various sources of 
ultraviolet (UV) radiation, like an active galactic nucleus (AGN) or a population 
of massive young stars (population III, see Carr, Bond and Arnett | 1984[ ) at 



the centre of the PG. So, the ultimate fate of the PGC clouds strongly depends 
- apart from other environmental conditions (more about this, later) - on the 
survival of H 2 . 

In the early phase of the PG an AGN is expected to form at its centre along 
with population III stars, because of the disruption of central PGC clouds. This 
indeed happens, for the cloud collision time is shorter than the cooling time in 
the central region of the PG. 

Thus, H 2 will be dissociated at galactocentric distances smaller than a cer- 



tain critical value R CT it- Following the analysis of Kang et al. [199C|, it is 



straightforward to evaluate R cr it- Consider first the case of a UV flux due to a 

6 Observe that — at variance with reactions (^) and ^ - reactions (^) do not require the 
presence of electrons and protons as a catalyst. 



central AGN. Then we find 



For typical luminosities up to Lagn — 10 45 erg s _1 , eq. (Q) yields R^+t* — 20 
kpc. On the other hand, when the UV dissociating flux is produced by massive 
young stars mainly located at the centre of the PG, the critical galactocentric 
distance turns out to be 

In eq. @ ~ 2 x 10 38 erg s _1 is the bolometric luminosity of a single BO 
V star. Assuming a total stellar luminosity L to t up to ~ 2 x 10 45 erg s _1 and 
a central number density uq up to ~ 10 3 kpc -3 , we find R* rit — 10 kpc. In 
conclusion, Hz should remain undissociated at galactocentric distances larger 
than 10 - 20 kpc. 

2.1 Globular clusters 

According to the foregoing analysis, in the inner galactic halo - that is for 
galactocentric distances smaller than 10 — 20 kpc - H% gets dissociated Q thus 
preventing any further cooling of the PGC clouds below T ~ 10 4 K. Therefore, 
these clouds remain for a long time in quasi-hydrostatic equilibrium at T ~ 10 4 
K (namely we have t coo i > tff during this period). In such a situation a 
characteristic mass scale gets imprinted on the PGC clouds by the gravitational 
instability, thereby resulting in a strongly peaked mass spectrum of the PGC 
clouds. In fact, Fall and Rees |1985| have shown that the Jeans mass of the 



7 T liis indeed occurs for a wide range of UV fluxes and PGC cloud densities (Kang et al. 
[1990|). 
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PGC clouds after the long permanence at T ~ 10 4 K is given (as a function of 
the galactocentric distance R) by ^] 

M PGC {R) ~ 5 x 10 5 (fi/kpc) 1 / 2 M , (9) 

while the PGC radius turns out to be 

r PG c(R) 20 (tf/kpc) 1 / 2 pc . (10) 

Observe that in the present case the propagation of sound waves erases all 
large-scale perturbations, leaving only those at small scales. 

Surely, this is not the end of the story, since the UV flux is expected to 
eventually decrease. So, after some time a nontrivial fraction of H2 should 
anyway form, causing a sudden drop of the PGC cloud temperature well below 
~ 10 4 K (the clouds now enter the regime t coo i < tff). What next happens is a 
rapid growth of the small-scale perturbations, which lead directly (in one step) 
- due to the thermal instability - to formation of stars inside the PGC clouds 



(Murray and Lin [ 1989 ) . It goes without saying that we expect supernova 
explosions to have occurred, thus giving rise to shocks inside globular clusters 
and, therefore, leading to a rather wide stellar IMF. In this way, globular clusters 
should form as they are observed today especially in the inner part of the galactic 
halo. Incidentally, the formation of the PGC clouds could have been delayed 
until after the Galaxy was enriched by metals due to population III stars, in 
this manner explaining the absence of globular clusters with primordial metal 
abundances and the radial gradient of metallicity in the galactic halo. 

8 We stress that the R-depende ncE in eqs. (0) and ( |lo[ ) holds for R < 20 kpc only (see 
the discussion in Vietri and Pesce ]1995| ), while for R > 20 kpc the mass of globular clusters 
tends rather to decrease. 
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2.2 Dark clusters of MACHOs and cold molecular clouds 

As we have seen, in the outer galactic halo - namely for galactocentric distances 
larger than 10 — 20 kpc - molecules are not dissociated due to the absence 
of a significant UV flux. Under the assumption of a quiet environment, this 
circumstance entails a very efficient cooling of the PGC clouds, whose state is 
therefore characterized by the condition t coo i < tff. Hence, the gravitational 
collapse is expected to occur as in the scenario of PSS. Specifically, both the 
temperature decrease and the density increase cause a substantial drop of the 
Jeans mass for the PGC clouds. As a consequence, they fragment into smaller 
and smaller clouds. This process stops when the clouds become optically thick 
to their own radiation, since then cooling manifestly gets uneffective. PSS have 
shown that such a situation occurs when the Jeans mass is as low as 10~ 2 M Q . 
We remark that the presence of rotation and magnetic fields in the PGC clouds 
would allow the Jeans mass to drop even further. Obviously, the fragmentation 
down to low-mass values is favoured if the initial gas has been metal enriched 
by population III stars (see Figure 7 in Palla and Stabler |1988| ). 



The result of the above picture is the formation of dark clusters containing 
compact objects with an IMF peaked in the range 10~ 2 — 10 _1 M Q . This value 
is close to the peak mass ~ 0.1 Af© of the present day IMF in the disk (Miller 
and Scalo |1979 ), but in our case we expect it to decrease much more rapidly 



for higher masses and thus to be more narrow around its mean value. So, the 
compact objects in question should be predominantly brown dwarfs, even though 
a fraction of M-dwarfs can be present. Actually, we want to stress that our 
considerations still hold true even if a large fraction of MACHOs consists of M- 



12 



dwarfs [] (as already mentioned, Kerins [1997] has shown that M-dwarfs clumped 
into clusters is a viable possibility to date). Notice that the expected lower limit 
of ~ 10~ 2 Mq is consistent with present bounds coming from microlensing data, 
as found by the EROS collaboration (Renault et al. 1997| ). 

What are the environmental conditions (besides H 2 survival) which permit 
the mechanism in question to work? First, turbulence should be negligible. In- 
deed, turbulence effects would make fragments collide and - as shown by PSS - 
this fact would increase the Jeans mass. Second, no sizable gravitational pertur- 
bations (shocks, supernova explosions, etc.) should be present for otherwise 
they would induce the gravitational collapse before cooling has succeeded in 
lowering the fragment Jeans mass down to the above-mentioned values. Be- 
cause these environmental conditions are likely to have occurred in the outer 
halo, the mass of the produced compact objects is expected to be close to the 
corresponding Jeans mass. 

Besides individual brown dwarfs, it seems quite natural to suppose that - 
much in the same way as it happens for ordinary stars - also in this case the 
fragmentation process should produce a substantial fraction of binary brown 
dwarfs. For reasons that will become clear later, they will be referred to as 
primordial binaries (we shall come back to this issue in the next Sections). It is 
important to keep in mind that the mass fraction of primordial binaries can be 



as large as 50% (Spitzer and Mathieu [|1980|). 

So, we are led to the conclusion that MACHOs consist of both individual 

9 The reader should keep this point in mind, in spite of our taking about brown dwarfs. Of 
course, quantitative changes in some results will occur, because M-dwarfs are more massive 
than brown dwarfs (this is especially true concerning Sects. 3 and 4). 

10 This circumstance is consistent with our assumption that the stellar IMF in the dark 
clusters is more sharply peaked than the IMF in the globular clusters and in the disk. 
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and binary brown dwarfs in this scenario (more about this, later). 

However, we don't expect the fragmentation process to be able to convert 
the whole gas in a PGC cloud into brown dwarfs. For instance, standard stellar 
formation mechanisms lead to an upper limit of at most 40% for the conversion 



efficiency (Scalo [ 1985 ). Therefore, a fairly large amount of gas - mostly H 2 
- should have been left over. What is its fate? At variance with the case 
of globular clusters, strong stellar winds are now manifestly absent, and so 
this gas should remain gravitationally bound in the dark clusters. Actually, 
this conclusion is further supported by the following arguments (provided dark 
clusters comprise a consistent fraction of the galactic dark matter). First, the gas 
cannot have diffused into the whole halo, for otherwise it would have been heated 
by the gravitational field to a virial temperature ~ 10 6 K, thereby becoming 
observable in the X-ray band - this option is ruled out by the available upper 
limits (Dickey and Lockman [ 1990| ) qj. Second, the alternative possibility that 
the gas wholly collapsed into the disk is also excluded, because the disk mass 
would then be of the order of the inferred dark halo mass. Now, the virial 
theorem entails that the temperature Tdg of a diffuse gas component inside a 
dark cluster is 

Accordingly, a large fraction of diffuse gas would presumably give rise to an 
unobserved radio emission. Thus, we conclude that the amount of virialized 
diffuse gas inside a dark cluster has to be low (it will henceforth be neglected). 
This circumstance implies in turn that most of the leftover gas should be in the 
form of self- gravitating clouds clumped into the dark clusters (since in this case 

11 Incidentally, the same argument also rules out a halo primarily made of unclustered brown 
dwarfs (as well as white and M-dwarfs). 
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the virial theorem applies to individual clouds). As we shall see later, there are 
good reasons to believe that the central temperature T m of the molecular clouds 
in question should be very low, in fact close to that of the CBR. Accordingly, 
the molecular cloud mass M m and median radius r m are related by the virial 
theorem as 

r ra ^4.8xl0- 2 ^ pc. (12) 

Presumably, the fraction of cluster dark matter in the form of molecular clouds 
should be a function of the galactocentric distance i?, depending on the en- 
vironmental conditions, like the UV flux and the collision rate for the PGC 
clouds. 

Before proceeding further, an important issue should be faced. Given the 
supposed existence of a large amount of gas in the dark clusters, one would 
expect a star formation process to be presently operative. However, things are 
not so simple. For - under the above environmental conditions - only stars with 
mass smaller than the cloud mass can be formed. Evidently, these stars are again 
either brown dwarfs or M-dwarfs. So, we see that undetected bright stars should 
not form in the dark clusters to the extent that our assumptions hold true. One 
might also wonder whether a sizable amount of gas is eventually left over. As 
already pointed out, we argue that this should be the case, for otherwise it 
would mean that the brown or M-dwarf formation mechanism should be much 
more efficient than any known star formation mechanism. Moreover, Gerhard 



and Silk [ 1996 1 have shown that the cluster gravitational field can stabilize the 
clouds against collapse. 

Unfortunately, the lack of any observational information about dark clusters 
would make any effort to understand their structure and dynamics hopeless, 
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were it not for some remarkable insights that our unified treatment of globular 
and dark clusters provides us. 

In the first place, it looks quite natural to assume that also dark clusters 
have a denser core surrounded by an extended spherical halo. For simplicity, we 
shall suppose throughout that the core density profile can be taken as constant. 
Moreover, it seems reasonable to imagine (at least tentatively) that dark clusters 
have the same average mass density as globular clusters. Hence, we obtain 

M \ 1/3 



rDC * °- 12 [-W ) pc ' (13) 

where Mdq and roc denote the mass and the median radius of a dark cluster, 
respectively. In addition, dark clusters - just like globular clusters - presumably 
stay for a long time in a quasi-stationary phase, with an average central density 
p*(0) slightly lower than 10 4 M Q pc -3 (which is the observed average central 
density for globular clusters). 

As a further implication of the present model, we stress that - at variance 
with the case of globular clusters - the mass spectrum of the dark clusters should 
be smooth, since the monotonic decrease of the PGC cloud temperature fails to 
single out any particular mass scale. As it will be shown in Sect. 3, dark clusters 
in the mass range 3 x 10 2 M Q < Moc < 10 6 M@ should survive all disruptive 
effects, and so we shall restrict our attention to such a mass range throughout. 

As far as dark clusters are concerned, we have seen that the brown dwarf 
mass is expected to lie in the range 10~ 2 — 10 _1 M Q . For definiteness (and with 
an eye to microlensing results) , we imagine that all individual brown dwarfs have 
the same mass to ~ 0.1 Mq. So, binary brown dwarfs are twice as heavy. As a 



consequence, the mass stratification instability (Spitzer [196E]) will drive them 



into the dark cluster cores, which then tend to be composed chiefly by binaries. 
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Furthermore, an average MACHO mass somewhat larger than ~ 0.1 M© can 
naturally be accounted for. 

Finally, let us consider molecular clouds. Since they also originate from the 
above-mentioned fragmentation process, we suppose (for definiteness) that they 
lie in the mass range 10~ 3 M Q <M m <10 _1 M . Correspondingly, eq. ( |l2| ) 
entails 4.8 x 10~ 5 pc<r m <4.8 x 10~ 3 pc and 2.7 x 10 10 cm~ 3 >n m >2.7 x 
10 6 cm -3 , respectively, where n m denotes the number density in the clouds. 

Before closing this Section, some comments are in order. There is little 
doubt that the foregoing considerations are qualitative in nature. Nevertheless, 
they provide nontrivial insights into several questions that arise in connection 
with the discovery of MACHOs. Specifically, a sharply peaked IMF in the range 
10~ 2 — 10 _1 Mq comes out naturally without having to invoke any new physical 
process. This fact explains the observed absence of a substantial amount of 
main-sequence stars inside the dark clusters. Therefore, the observed absence 



of a large number of planet-like objects in the halo (Renault et al. [1997]) is 
automatically explained. Furthermore, we can understand why brown dwarfs 
clumped into dark clusters copiously form in the outer halo, but neither in the 
inner halo nor in the disk. Indeed, the different stellar content of these regions is 
here traced back to the different environments in which the same star formation 
mechanism operates. It goes without saying that various issues addressed above 
require further investigations. 

3 Dynamical constraints on dark clusters 

As we have seen, MACHOs are clumped into dark clusters when they form in 
the outer galactic halo. Still, the further fate of these clusters is quite unclear. 
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For, they might either evaporate or drift towards the galactic centre. In the lat- 
ter case, encounters with globular clusters might have dramatic observational 
consequences and dynamical friction could drive too many MACHOs into the 
galactic bulge. So, even if dark clusters are unseen, nontrivial constraints on 
their characteristic parameters arise from the observed properties of our galaxy. 
Moreover, in order to play any role as a candidate for dark matter, MACHOs 
must have survived until the present in the outer part of the galactic halo. Fi- 
nally, it is important to know whether MACHOs are still clumped into clusters 
today, especially because an improvement in the statistics of microlensing ob- 



servations permits to test this possibility (Maoz [ 1994 1 , Metcalf and Silk [1996]). 
We remark that previous work on dynamical constraints on clusters of brown 



dwarfs (Carr and Lacey |1987|, Carr [1994|, Moore and Silk |199S], Gerhard and 



Silk [1996 1, Kerins 1996 ) rests upon the hypothesis of an initial dark cluster 
distribution that extends inside the inner part of the Galaxy all the way down 
to the centre. Hence, a novel ab initio analysis is required within the present 
scenario. 

We begin with the usual assumption that the halo dark matter density is 
modelled by an isothermal sphere with density profile 



p{R) = pa 



R 2 



(14) 



where a ~ 5 kpc is the halo core radius, Rq ~ 8.5 kpc is our galactocentric 
distance and po ~ 6.5 x 10~ 25 g cm~ 3 denotes the local dark matter density 
corresponding to the one-dimensional velocity dispersion a ~ 155 km s _1 of 

12 As mentioned in Sect. 1, this model may not provide the best description of the Galaxy. 
However, it will still be used in the present Section since otherwise the ensuing discussion 
would become exceedingly complicated. Furthermore, it should be kept in mind that the 
various uncertainties affecting the dark cluster properties anyway make the following results 
reliable only as order-of-magnitude estimates. 
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the halo. Furthermore, we shall suppose for definiteness that the age of the 
Universe is to ~ 10 10 yr. Our treatment of encounters rests upon the diffusion 
approximation and we will follow rather closely the analysis of Binney and 



Tremaine [1987 



3.1 Dynamical friction 

Dark clusters are subject to dynamical friction as they orbit through the Galaxy, 
which makes them loose energy and therefore spiral in toward the galactic centre. 
Assuming (for illustrative purposes) that a dark cluster moves with velocity v on 
a circular orbit of radius R, the drag brought about by the background density 
p(R) is given by 

F{R) = - 0.43 22 p(R) 15 

v 

where In A is the usual Coulomb logarithm, whose value in the present case is 
In A ~ In (fg^) - 24.3 - ln(M DC /M G ) , (16) 



with R typ ~ 20 kpc, v ~ \/2a. Using eq. (|l4|), eq. ([15]) becomes 

F(R) ~ - 5.1 x 10 14 lnA g cm s" 2 , (17) 

where i?5 is in units of 5 kpc. Accordingly, the equations of motion entail that a 
dark cluster originally at galactocentric distance R will be closer to the galactic 
centre today by the amount 

AR(R) * 1.85 X «r» ("SSL) »*-HM*ol*>) kp , (18) 
V m & j r 5 + r 5 

Keeping in mind that in our model R > 10—20 kpc and Mdc < 10 6 M© 
(see later), we see that Ai?<5.8 x 10" 2 kpc. Therefore, dark clusters are still 
confined in the outer galactic halo. As a consequence, encounters between dark 
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and globular clusters as well as disk and bulge shocking of dark clusters are 
dynamically irrelevant, as long as they move on not too highly elongated orbits 
(in this way an effective circularization of the orbits is achieved). 

3.2 Encounters between dark clusters 

Encounters between dark clusters may - under the circumstances to be analyzed 
below - lead to their disruption. As an orientation, we notice that an estimate 
of the one-dimensional velocity dispersion a* of MACHOs and molecular clouds 
within a dark cluster is supplied by the virial theorem and reads 

1/3 

a. ~ 6.9 x 10- 2 ( km s" 1 , (19) 

where eq. ([l3|) has been used. Because in the present scenario M DC <10 6 M Q 
(see later), we get cr» <6.9 km s _1 . Therefore, the one-dimensional velocity 
dispersion of dark clusters - which we naturally suppose to be just a - is much 
larger than ov Hence it makes sense to work within the impulse approximation, 
whose range of validity is established more precisely by the condition Mjjc *C 
10 10 Mg, which is evidently always met in our model. 

In order to proceed further, we denote by AE the change of the internal 
energy of a dark cluster in a single encounter. Then (following Binney and 
Tremaine 1987| |) we find that encounters with impact parameter b in the range 



&min < b < & max increase the cluster's energy at the rate 



E{R) ~ ^ ' / dvv A e~ v / 4ff / dbbAE (20) 

'0 Jfcmin 



where v and Udc{R) are the cluster velocity and number density (in the halo), 
respectively. We let 7 stand for the fraction of halo dark matter in the form of 
dark clusters, so that we have noc(R) = 1p{R)/Mdc, with p(R) given by eq. 
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(p"4|). Accordingly, eq. ( p0| ) becomes 

- 7"7e TTTTT" 7777-7 / ^ 3 e^ 2/4CT2 jf d& 6 AS . (21) 



20F ctGM D c R 2 + a 2 J 



Now, a natural definition of the time required by encounters to dissolve a cluster 
is provided by td{R) = E\,- m &l E(R), where the binding energy i?bind is expressed 
in terms of the cluster properties as -Ebind — 0.2 GMj^q/tdc- 

3.2.1 Distant encounters 

Let us consider first distant encounters. Correspondingly, AE is to be evaluated 



in the tidal approximation (Spitzer [ 1987 1 ) and reads presently 



AE * - 2?f DC ■ (22) 

We insert this quantity into eq. (pl|), neglecting the term (& m i n /6 ma x) 2 in the 
ensuing expression. Experience with a similar treatment for globular clusters 
suggests to choose & m ; n ~ rpc- Correspondingly, we find 

4.7 x 10 11 / M@ \ 1/3 



^-^"f-^J ^ + Dyr, (23) 
where use of eq. ( |l3| ) has been made. Assuming R > 10 — 20 kpc and keeping 
in mind that 7 < 1 and Moc < 10 6 Mq (see later), we get that for all the dark 
clusters in question td(R) exceeds the age of the Universe. 

3.2.2 Close encounters 

In order to deal with close encounters, dark clusters have to be regarded as 
extended objects. As in the case of globular clusters, this task is most simply 
accomplished by modelling the dark clusters by means of a Plummer potential 
with core radius a. Correspondingly, AE is found to be 
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As before, we insert this quantity into eq. (pl|), assuming now 6 m i n ~ and 
6 max ~ Tdc- ln this way, we obtain 

. . 2.4 x 10 11 / a \ 2 /pc\, , 

t d {R ~ M- (i? + 1 yr. (25) 

7 \tdcJ \rDC J 

We proceed by recalling that 

/ 3M DC x 1/3 



V47rp >t (0) 



(26) 



where p*(r) denotes the dark cluster mass density. Taking p,(0) ~ 10 4 M 
pc~ 3 (as for globular clusters today) and using eq. (|T^), we can rewrite eq. 



I.U" f Mq \ 1/3 2 



as 

t d {R) — ~~ — ( ~ip^~ ) (^5 + 1) yr. (27) 

7 \ Ad D c J 

Assuming R > 10 — 20 kpc and remembering that 7 < 1, we are led to the 
conclusion that dark clusters are not disrupted by close encounters provided 

M DC <10 6 Mq. □ 

3.3 Evaporation 

Various dynamical effects conspire to make dark clusters evaporate within a 
finite time. Relaxation via gravitational two-body encounters leads to the escape 
of MACHOs approaching the unbound tail of the cluster velocity distribution. 
Tidal truncation due to the galactic gravitational field enhances this process. A 
more substantial effect is caused by the gravothermal instability, when the inner 
part of the dark clusters contracts (core collapse) and the envelope expands. 

Below, we shall address these issues separately. As is well known, a key-role 
in such an analysis is played by the relaxation time 

^ I °- 34 G ^/;in(0.4A0 ^ 



13 It should hardly come as a surprise that close encounters yield a more stringent bound on 
Mjjc than distant encounters. 
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where N is the number of MACHOs per cluster. As in the case of globular 
clusters, p*(r) is expected to vary by various orders of magnitude in different 
regions of a single dark cluster - this dependence obviously shows up in i re iax(?~)- 
Therefore, for reference purposes, it is often more convenient to characterize a 
dark cluster by a single value of the relaxation time. This goal is achieved by 



introducing the median relaxation time (Spitzer and Hart |1971|) 



6.5 x 10* ( M DC \ 1/2 ( Mq \ fr DC ^ 3/2 



tih= Hom\^M^) {—){^) yr - (29) 

Explicitly, using eq. ( |l5| ) together with eq. ( |T^ ) and m ~ 0.1 M , eq. ( p8| ) 
takes the form [^] 

ireiax r ~ 5 x 10 — — -—— — — yr . 30 

V M Q J \ p*(r) J 1.4 + ln(M DC /M ) 

In the same fashion, eq. (p9| ) becomes 

,5 ( M DC 



t ^ 8Xl ° 5 ^J » tM « K /M 8 ) -' < 31 » 
3.3.1 Spontaneous evaporation 

As is well known, any stellar association evaporates spontaneously q] within a 
finite time owing to relaxation via gravitational two-body encounters. Specifi- 
cally, a single close encounter between two MACHOs can leave one of them with 
a speed larger than the local escape velocity. So, the MACHO under consider- 
ation gets ejected from the dark cluster. We find for the ejection time (Hcnon 



[1969 ) 



i ej ~ 1.1 x 10 3 ln(0.4A) t rh ~ 9 x 10 8 ( ) yr. (32) 



M 



14 For simplicity, we neglect here the fact that binaries have a mass larger than individual 
brown dwarfs. Furthermore, given the logarithmic JV-dependence, we can safely take N ~ 
M D c/m ~ 10 (M DC /M e ). 

15 We like to call in this way the evaporation process which is neither induced by external 
perturbations nor specific of the gravitational interactions. 
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Alternatively, several more distant weaker encounters can gradually increase the 
energy of a given MACHO, until a further weak encounter is sufficient to make 
it escape from the cluster. In this case, the evaporation time turns out to be 



(Spitzer and Thuan 1972 ) 



W P * 300 t, h * 2.4 x 10- (*g>) l4 + m ; wMq) (33) 

Since i e j is anyway longer than £ G vapj w ^ shall focus our attention on the latter 
quantity. Hence, by demanding that i CV ap should exceed the age of the Universe 
we conclude that dark clusters with Mnc > 3 x 10 2 Mq are not yet evaporated. 

3.3.2 Tidal Perturbations 

Dark clusters - just like globular clusters - are tidally disrupted by the galactic 
gravitational field unless roc is smaller than their tidal radius. So, the survival 
condition reads 

where R should be understood here as the perigalactic distance of the dark 
cluster and Mg{R) denotes the Galaxy mass inside R. From eq. ( |l4| ) we obtain 

M G (R) ~ 5.5 x 10 10 R 5 (l - R^ 1 arctani? 5 ) M Q , (35) 



and - on account of eq. (13) - eq. (EM) becomes 



R 5 [l - i?«r 1 arctani? 5 ]~ 1/2 > 0.047 , (36) 

which is always satisfied for R > 10 — 20 kpc. Thus, the dark clusters under 
consideration are not tidally disrupted by the galactic gravitational field. 
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3.3.3 Core collapse 

An important role in the considerations to follow is played by core collapse. 
It is by now well-established that the initial stage of this process is triggered 
by evaporation, which leads to the shrinking of the core as a consequence of 
energy conservation. Numerical Fokker-Planck studies of the early phase of 
core collapse have shown that the dynamics of the cluster is correctly described 



by a sequence of King models (Cohn 1980 ). However, once the cluster density 



reaches a certain critical value, core collapse gets dramatically accelerated by 



the gravothermal instability (Antonov [ 1962 1 , Lynden-Bell and Wood [1968]). 



Indeed, the negative specific heat of the core implies that the internal velocity 
dispersion tr* increases - thereby enhancing evaporation - as the average kinetic 
energy decreases by evaporation itself. Moreover, the unbalanced gravitational 
energy makes the core contract, and so its density rises by several orders of 
magnitude in a runaway manner. Numerical simulations show that the central 
velocity dispersion and the number of stars in the core iV* scale as 

<r. ~ p»(0)°-° 5 , (37) 

N* ~ /9*(0)-°- 36 , (38) 



respectively (Binney and Tremainc [1987]). Incidentally, the somewhat surpris 



ing slow rise of a* in eq. (^) is due to the large mass-loss from the core, as it 
follows from eq. (pq). 

When does the gravothermal instability show up? Unfortunately, a clear-cut 
answer does not exist, since the corresponding time tai depends on how clusters 



are modelled as well as on their concentration (Quinlan 1 1996 ] ) . Manifestly, the 



lack of observational data about dark clusters makes a precise determination of 
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tGi impossible. So, the best we can do is to suppose that dark clusters behave 
just like globular clusters as far as core collapse is concerned. In this way, we are 



lead to the order-of-magnitude estimate (Binney and Tremaine [1987]) (more 
about this, later) 

'""^"'■''(t) 1.4 + 1 „(M DC /M e ) yr - (39) 
Comparing now tQj with the age of the Universe, we conclude that the dark 
clusters with Mdc <5 x 10 Mq are expected to have started core collapse. 

As we said, the central density grows dramatically during the second stage 
of core collapse, and so the central relaxation time gets shorter and shorter. 
Detailed studies of the gravothermal instability have shown that - if nothing 
opposes to collapse - the time needed to complete core collapse t co ii starting from 
an arbitrary time t goes like i re iax(0), with the latter quantity computed for the 
particular value taken by p*(0) at t. More specifically, computer simulations of 
globular cluster dynamics entail t co \\ — 330 trelax(O) (Cohn jjggjfl ). 



Because of the huge increase in the central density, close two-body encounters 
lead to the formation of bound binary systems by converting enough kinetic 
energy into internal energy (tidal capture). As we will see in Section 4, binary 
brown dwarfs are produced in this way in the dark cluster cores during the 
early phase of core collapse. These binaries - which will be referred to as tidally- 
captured binaries and happen to be very hard - play a crucial role in this context, 
since they ultimately stop and reverse the collapse. Schematically, the argument 
goes as follows. Owing to the fact that hard binaries necessarily get harder 



in collisions with individual stars (Heggie [1975Q, the internal binding energy 
released by a binary is transformed into kinetic energy of both the star and 
the binary. Actually, the exchanged energy is so large that they both leave 

2G 



the cluster. However - at variance with evaporation - the kinetic energy (per 
unit mass) of the cluster is now unaffected, whereas mass ejection obviously 
increases the potential energy. That is, the binding energy given up by the 
binaries ultimately becomes gravitational energy of the core. As a result of 
the unbalanced kinetic energy, the core starts expanding. Moreover, because 
of the negative specific heat the increased potential energy makes cr* decrease, 
thereby slowing down mass ejection. In this manner, core collapse gets halted 



and reversed (Spitzer 1987 ) 



As a matter of fact, the presence of binaries in appreciable amount can also 
modify to some extent the standard scenario of core collapse as outlined above. 
Indeed, numerical Fokker-Planck simulations have shown that in this case the 
collapse is also driven by the mass stratification instability. As a consequence, 
the collapse proceeds faster and starts before than in the few-binary case (the 
latter point makes eq. (|3^) more plausible than it might appear at first sight). 
This phenomenon is found to occur both for tidally-captured binaries (Statler, 



Ostriker and Cohn [1987]) and for primordial binaries (Spitzer and Mathieu 



Hi 9801 ) 



3.4 Discussion 

What the above analysis shows is that dark clusters within the mass range 
3 x 10 2 Mq < Mdc < 10 6 Mq should have survived all disruptive effects aris- 
ing from gravitational perturbations and are nowadays expected to populate the 
outer part of the galactic halo. In addition, the clusters with 3x 10 2 M Q < Mdc < 5x 
10 4 Mq should undergo core collapse. Unfortunately, it is practically impossible 
to predict the further fate of those dark clusters which are in the post-collapse 



27 



phase today. A priori it seems natural to imagine that bounce and subsequent 
reexpansion should follow core collapse (Cohn and Hut [ 1984 1 , Heggie and Ra- 



mamani | 1989fl ). And perhaps also a whole series of core contractions and ex- 



pansions can take place, giving rise to the so-called gravothermal oscillations 
(Bettwieser and Sugimoto 1984| ). However, this conclusion crucially depends 



on the unknown model which correctly describes dark clusters. For instance, in 
the case of tidally-truncated models the cluster gets completely destroyed within 
a finite time (Stodolkiewicz | 1985|| , Ostriker, Statler and Lee [ 1985 1 ) . Moreover, 



what certainly happens in either case is that the number of MACHOs in the core 
monotonically decreases with time. So, an unclustered MACHO population is 
expected to coexist with dark clusters in the outer galactic halo (unless all dark 
clusters have M DC >5 x 10 4 M Q ) - detection of unclustered MACHOs would 
therefore not rule out the present scenario. 

4 MACHOs as binary brown dwarfs 

As already pointed out, it seems natural to suppose that a fraction of primordial 
binary brown dwarfs - possibly as large as 50% in mass - should form along with 
individual brown dwarfs as a result of the fragmentation process of the PGC 
clouds. Subsequently, because of the mass stratification instability primordial 
binaries will concentrate inside the dark cluster cores, which are therefore ex- 
pected to be chiefly composed by binaries and molecular clouds. In addition - 
as far as dark clusters with M^c < 5 x 10 4 M Q are concerned - a population of 
tidally-captured binary brown dwarfs ought to form in the dark cluster cores 
owing to the increased central density caused by core collapse. Thus, a large 
fraction of binaries should be present inside the dark cluster cores at a late 
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stage of their evolution. Below, we will try to make the discussion of this issue 
as much quantitative as possible. 

4.1 Survival and hardness of binary brown dwarfs 

The first question that has to be addressed is whether a binary brown dwarf - 
produced by whatever mechanism long ago - survives up until the present. To 
this end, we recall that a binary system is hard when its binding energy exceeds 
the kinetic energy of field stars (otherwise it is soft). In the present case, binary 
brown dwarfs happen to be hard when their orbital radius a obeys the following 
constraint 

/ M \ 2 / 3 

a<1.4xl0 12 — 2-1 km. (40) 



\M DC/ 

As is well known, soft binaries always get softer whereas hard binaries always 



get harder because of encounters with individual stars (Heggie [ 1975 1 ) . So - were 
individual-binary encounters the only relevant process - we would conclude that 
hard binary brown dwarfs should indeed survive. However, also binary-binary 
encounters play an important role in the dark cluster cores, where binaries are 
expected to be far more abundant than individual brown dwarfs. Now, in the 
latter process one of the two binaries gets often disrupted (this cannot happen 
for both binaries - given that they are hard - while fly bys are rather infrequent), 
thereby leading to the depletion of the binary population. We will address this 
effect in Sect. 4.3, where we shall find that for realistic values of the dark cluster 
parameters the binary break up does not take place. 

4.2 Tidally-captured binary brown dwarfs 

As far as globular clusters are concerned, it is nowadays well known that the 
most efficient mechanism for late binary formation is dissipative tidal capture in 
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the core (Fabian, Pringle and Rees [ 1975 1 , Press and Teukolsky 1977 , Lee and 
Ostriker | 1986| ). Hence, we expect a similar situation to occur for dark clusters 
as well. 

Let us now analyze this phenomenon in a quantitative fashion. As a first 
step, we observe that the radius r* of a brown dwarf of mass ~ O.1M is 
r„ ~ 0.7 x 10 5 km (Saumon et al. | 1996| ). Furthermore, the Safronov number 



(Binney and Tremaine [ 1987 1 ) is 

e= Gm 



2 ol r* 



2 x 10' 



( Me \ 



2/3 



(41) 



\M DC J 

which turns out to be anyway much larger than one. Within this approximation, 
the time for brown dwarf tidal capture can be written as (Lee and Ostriker 



[ |1986| ) 



tad - 10 



12 / 10 5 PC~ 3 



12 /i? V" fM t 



0.9 / , r \ 1.1 





yr, 



(42) 



niBD{0) J V100 km s" 1 , 

where 71jbd(0) is the number density of individual brown dwarfs in the core P*| 
Obviously, we have n/B_o(0) ~ fiBD p*(0)/m, with fiBD denoting the mass 
fraction of individual brown dwarfs in the core. From eq. (|37| ) we see that cr* 
increases very slightly, and so core collapse effects on er* can safely be neglected. 
Accordingly, we can rewrite eq. (E2h as 



0.4 



/ ^ 1.6 x 10 14 / M Q pc~ 3 \ ( M D c 



fiBD V P*(0) J \ M e J ^ (43) 

where eq. (19h has been used. Comparing now i t id with the age of the Uni 



verse, we get that practically all individual brown dwarfs in the core are tidally 

16 Although eq. (|42|) has been derived for main-sequence stars, it looks plausible that it 
applies also t o brow n dwarfs, since it is not very sensitive to the particular stellar model (Lee 
and Ostriker [h.986]). 
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captured into binaries provided ^| 



,(«)>^f¥l (44) 



flBD V M Q 

According to the above assumptions, we expect p*(0) ~ 10 4 M pc~ 3 just 
before core collapse. Therefore, we see that tidal capture requires an increase 
of p*(Q) by a factor in the range 31/ Jibd ~ 242/ Jibd, corresponding to Mdg 
in the range 3 x 10 2 M — 5 x 10 4 Mq. Thus, the formation of tidally-captured 
binaries would occur during the early phase of core collapse (the same conclusion 
was reached in a different way by Statler, Ostriker and Cohn 1987| for globular 
clusters). However, this conclusion depends on the fractional (mass) abundance 
fpB of primordial binaries in the core, since /ibd necessarily gets small for 
large f PB - 

Next, we compute the (average) orbital radius of tidally-captured binary 
brown dwarfs following the procedure outlined by Statler, Ostriker and Cohn 



[1987|. Correspondingly, we find a ~ 2.5 x 10 5 km (this value is practically 
independent of M^c)- As a consequence, we see that they are so hard that 
condition ( |40| ) is always abundantly met. 

Let us finally try to estimate the fractional abundance of tidally-captured 
binary brown dwarfs in the dark cluster cores soon after their formation, namely 



when inequality (44) just starts to be satisfied (their total number at this stage 
will be denoted by Ntcb)- Thanks to eq. ( |38l ) we easily get 

WTCB ,,3^(^) 08, (^), ,45, 



17 More precisely, we are demanding that the rate for tidal capture times the age of the 
Universe should exceed one. Notice that the former quantity is one-half of Cl, and so we 
require that 2t t i c j should be sm aller than 10 10 yr (we are actually following the same procedure 
used by Press and Teukolsky [1977 ). 
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where M c denotes the core mass just before core collapse (corresponding to 
p*(0) ~ 10 4 Mq pc~ 3 ). Denoting further by N\°g D the total number of individ- 
ual brown dwarfs in a dark cluster before core collapse, we find 



Ntcb _ „ on f0 .36 

/ m q v v M c 



Realistically, even in the extreme case of a fully baryonic halo we expect Jibd < 0.3 
(of course, a large fps would imply a fiBD considerably smaller than that). 
Moreover, M c /Moc sensitively depends on how the core is defined, but the 
analogy with globular clusters entails that it should anyway be less than 20%. 
Accordingly, we conclude that the fractional abundance of tidally-captured bi- 
nary brown dwarfs should not exceed 1% (again in remarkable agreement with 



the result of Statler, Ostriker and Cohn [1987] for globular clusters), and so they 



are irrelevant from the observational point of view. 
4.3 Primordial binary brown dwarfs 

Primordial binaries are a very different story. For, not only are they expected to 
be much more abundant than tidally-captured binaries, but in addition they are 
presumably much less hard, since all values for their orbital radius consistent 
with condition ( ffOj ) are in principle to be contemplated. So, hardening effects 
ought to play a crucial role for primordial binaries (as it can be guessed on 
intuitive grounds, effects of this kind turn out to be totally negligible for tidally- 
captured binaries - this will be demonstrated later on). 

4.3.1 Collisional hardening 

Let us begin by considering collisional hardening, namely the process whereby 
hard binaries get harder in encounters with individual brown dwarfs. We recall 
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that the associated average hardening rate (Spitzer & Mathieu 1980) reads 
presently 

E * - 2.8 ° 2m3 UlBDi0) . (47) 

(7* 



On account of the definition of rz.jszj(0) and of eq. fll9|), eq. ( |47| ) becomes 

Observe that a characteristic feature of collisional hardening is that is inde- 
pendent of hardness, and so it is time-independent. 

As is well known, the internal energy of a binary is E = —Gin 2 /2a, which 
yields 

■ Gm 2 . , . 

E=^a. (49) 

Hence, by combining eqs. (E§) and ( |49| ) together and integrating the ensuing 
expression we get 



flBD fej Vm^^) [v) ' 



^^ + 1.3xl^ u k,(-! (^— ")(-) ■ 

where ai stands for the initial orbital radius, whereas a 2 denotes the orbital 
radius after a time £21- 

Assuming momentarily that no other hardening mechanism is operative and 
taking t 2 i equal to the age of the Universe, we find that the present orbital 
radius of a binary brown dwarf is given by 



^^ + 1 .3xlO-»/« I ,(-^) V3 ( 1 ^_ 
0,2 ai \M D c J \M Q pc 6 



Of course, collisional hardening works to the extent that a 2 becomes consider- 
ably smaller than a\. Correspondingly, we see from eq. ( ]5l| ) that this is indeed 
the case provided 

,-x {Mdc\ 1/3 (M q pc- 3 



-issxitfva.l-jjfj {-^r' km - ,52 » 
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Thanks to eq. (]5jj), eq. ( |52| ) yields in turn 

fl 2 - 8 x 10 / Jfl23 — — - km . (53) 



\M Q J V P.(0) 
Physically, the emerging picture is as follows. Only those binaries whose 
initial orbital radius obeys condition ( |52| ) undergo collisional hardening, and 
their present orbital radius turns out to be almost independent of the initial 
value. We can make the present discussion somewhat more specific by noticing 
that our assumptions strongly suggest p*(0) < 1O 4 M pc -3 , in which case both 
cq. ( j52|) and eq. (|33j) acquire the form 

1/3 

"i^SxlO 5 /^^) km. (54) 

Evidently, very hard primordial binaries - which fail to meet condition ( |54"|) 

- do not suffer collisional hardening, and the same is true for tidally-captured 
binaries. 

4.3.2 Prictional hardening 

As we are going to show, the presence of molecular clouds in the dark cluster 
cores - which is indeed the most characteristic feature of the model in question 

- provides a novel hardening mechanism for binary brown dwarfs. Basically, 
this is brought about by dynamical friction on molecular clouds. 

It is not difficult to extend the standard treatment of dynamical friction 



(Binney & Tremaine [1987]) to the relative motion of the brown dwarfs in a 
binary system which moves inside a molecular cloud. For simplicity, we assume 
that molecular clouds have a constant density profile p m . In the case of a circular 
orbit F 8 ], the equations of motion imply that the time 4i needed to reduce the 



18 Indeed, the circularization of the orbit is achieved by tidal effects after a few periastron 
passages (Zahn 1987). 
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orbital radius a of a binary which moves all the time inside molecular clouds 
from a\ down to ai is 



4? ~ 0.17 



( m \ 1/2 1 / -3/2 -3/2x 

Kg) ^TnA (fl2 - fll } ' (55) 

where the Coulomb logarithm reads 

In A ~ \n(r m v 2 c /Gm) ~ ln(r m /ai) , (56) 

with i> c denoting the circular velocity (approximately given by Kepler's third 
law). Manifestly, the diffusion approximation - upon which the present treat- 
ment is based - requires that the orbital radius of a binary should always be 
smaller than the median radius of a cloud. As we are concerned henceforth with 
hard binaries, a\ has to obey condition flip]). On the other hand, a% is the larger 



value for the orbital radius in eq. (55). So, we shall take for defmiteness - in 



the Coulomb logarithm only - a% ~ 1.4 x 10 12 (Mq/Mdc) 2 ^ 3 km. In addition, 
from eq. ( |l2] ) we have 

p m ~ 2.5 f ^ M pc~ 3 . (57) 



Hence, putting everything together we obtain 

^ a/2 ,« 3/2 + 2xl „-»H-.(-)7«), ,58, 
a-2 J \ ai J \r m J \ yr J 

having set 

S = [3 + ln(r m /pc) + O.71n(M flC /M )]- 1 . (59) 
Specifically, the diffusion approximation demands S > 0, which yields in turn 

r m > 5 x 10- ^L) pc . (60) 

Observe that for Mdc <2.1 x 10 4 M Q this constraint restricts the range of 
allowed values of r m as stated in Sect. 2.2. 
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Were the dark clusters completely filled by clouds, eq. (58) would be the 
final result. However, the distribution of the clouds is lumpy. So, if we want to 
know the orbital radius ai after a time ti\ we have to proceed as follows. First, 
we should compute the fraction of the time interval in question t2i spent by 
a binary inside the clouds. Next, we have to re-express t$ in eq. (^8|) in terms 
of ^2i- This goal will be achieved by the procedure outlined below. 

Keeping in mind that both the clouds and the binaries have average velocity 
v ~ -y/Scr* (for simplicity, we neglect the equipartition of kinetic energy of the 
binaries) it follows that the time needed by a binary to cross a single cloud is 

As an indication, we notice that for r m ~ 10 -3 pc (M m ~ 2 x 10~ 2 Mq) and 
Mdc — 10 5 Mq we find t m ~ 1.2 x 10 2 yr. Therefore, frictional hardening 
involves many clouds. Specifically, during the time £21 the number of clouds 
crossed by a binary is evidently 



N„ 



Let us now ask how many crossings of the core are necessary for a binary to 
traverse N m clouds. To this end, we proceed to estimate the number of clouds 
N c encountered during one crossing of the core. Describing the dark clusters by 
a King model, we can identify the core radius with the King radius. Evidently, 
the cross-section for binary-cloud encounters is Trr^ and so we have 

N ^{^Mm)'* ncL{QUrl > (63) 

with nci(O) denoting the cloud number density in the core. Thanks to eq. (pjj), 
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we can write 

ncUO) - foz.^ - 4.8 x 10-Vc, ( ^) ( jf^) pc" 3 , (64) 

where fcL denotes the fraction of core dark matter in the form of molecular 
clouds. Correspondingly, eq. ( |63|) becomes 

on account of eq. (|l9|). So, the total number of core crossings N cc that a binary 
has to make in order to traverse N m clouds is 

^(^fte) 2 (f)- « 

Because the core crossing time is 

it follows that the process under consideration takes a total time N cc t cc . How- 
ever, this time is by definition just t%]_. So, we get 

t 21 ^ cc ^,9.8 / -(^^)(^) 2 4°), (68) 

which is the desired relationship between t^x and t^x • We are now in position 
to re-express eq. (p8|) in terms of t2x- Accordingly, we get 

In order to quantify the effect of frictional hardening, we may proceed much 
in the same way as in the case of collisional hardening, neglecting however the 
latter effect for the moment. Specifically, taking t^x in eq. ([39]) equal to the age 
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of the Universe, we find that the present orbital radius of a binary brown dwarf 
is given by 

Manifestly, frictional hardening is operative to the extent that a-i becomes con- 
siderably smaller than a\. Accordingly, from eq. (|7^) we see that this is indeed 
the case provided 

kn - (71) 

Owing to eq. (|7l|), eq. ( |70| ) entails in turn 

O2 ,l,, 1 0../ c -r^(^) 2/3 to. (72, 

Physically, only those binaries whose initial orbital radius satisfies condition 
( frlj ) are affected by frictional hardening, and their present orbital radius turns 
out to be almost independent of the initial value. We can make the present 
discussion somewhat more specific by noticing that our assumptions strongly 
suggest /9*(0) < fO 4 M Q pc~ 3 , in which case both eq. (fn]) and eq. ( |72] ) acquire 
the form 

ai, 2 > 2.8 x 10 s f c l /3 Z 2/3 km . (73) 

Evidently, very hard primordial binaries - which violate condition (|73| ) - do not 
suffer frictional hardening, and the same is true for tidally-captured binaries. 

4.3.3 Present orbital radius of primordial binaries 

Which of the two hardening mechanisms under consideration is more effective? 
A straightforward implication of eqs. (^) and (|7^) is that frictional hardening 
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is more efficient than collisional hardening whenever it so happens that 

Unfortunately, the occurrence of various dark cluster parameters in condition 
( [74]) does not permit a sharp conclusion to be drawn, but in the illustrative 
case Mdc ^ 10 5 M o ,r m ~ 10~ 3 pc, f C L ^ 0.5 and p*(0) ~ 1O 3 M pc~ 3 we 
get fiBD < 0.5. As we expect in the cores f IB D < Ipb and f CL ~ f PB , we 
see that frictional hardening plays the dominant role. Moreover, from eqs. ( |53| ) 
and (|7^) it also follows that the effectiveness of collisional hardening decreases 
for smaller values of /j bd ■ Finally, the fairly slow dependence of condition (Q) 
on Mdc and p*(0) makes our conclusion rather robust. 

Having shown that collisional hardening can effectively be disregarded, we 
get that the present orbital radius of primordial binary brown dwarfs is actually 
given by eq. (|72|), as long as condition ( |7l| ) is met. Taking again the above 
particular case as an illustration, we find 02 — 7 x 10 s km, which is of the 
same order as the Einstein radius for microlensing towards the LMC (Gaudi and 



Gould [ 1997 1 ) . Therefore, we argue that not too hard primordial binaries can be 
resolved in future microlensing experiments with a more accurate photometric 
observation, the signature being small deviations from standard microlensing 



light curves (Dominik |1996| ) 



4.3.4 Gravitational encounters 

We are now in position to take up the question concerning the survival of bi- 
nary brown dwarfs against gravitational encounters. As already pointed out, 
individual-binary encounters are harmless in this respect - since hard binaries 
are considered throughout - and so we shall restrict our attention to binary- 
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binary encounters. 

As a first step, we recall that their average rate T (Spitzer & Mathieu 1980) 
can presently be written as 

„ Gma 

r ~ a , (75) 

(7* 

with a ~ 13. Correspondingly, on account of eq. ( [l9| ) the reaction time in the 
dark cluster cores turns out to be 

with /3 ~ 6.7//ps. 

Now, if no hardening were to occur - which means that the orbital radius 
a would stay constant - the binary survival condition would simply follow by 
demanding that t reac t should exceed the age of the Universe. However, hard- 
ening makes a decrease, and so t reac t increases with time. This effect can be 
taken into account by considering the averge value < t react > of the reaction 
time over the time interval in question (to be denoted by T), namely 

1 f T 

^ treact ^ — 7p / ^ ^react • (^7) 
1 JO 

In order to compute < t react >, the temporal dependence a(t) of the binary 
orbital radius is needed. Because frictional hardening plays the major role, 
the latter quantity is evidently supplied by eq. ([39]) ■ Setting for notational 
convenience t = i2i and a(t) = a-i-, eq. (|69|) yields 



a{t) 



km V V2 , 2.1 xlO" 27 -" 1 /^ PM U * 



-2/3 



km . (78) 



Combining eqs. (|7^ ) and ( |78| ) together and inserting the ensuing expression into 
eq. (0), we get 

/ ji r \ 1/3 /l \ 5/2 / „, r -3\ 2 



<t reQCt > ~ 1.9 x 10 46 f P lf c l S (^^J 



T / \ ai J \ p*(0) 
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1 + 2.1 x W-^E- 1 fcL 



P*(0) 



T 



3/2 



5/3 



v Mqpc 3 J \yr J Vkm7 
Let us now require < t react > to exceed the age of the Universe (taking evidently 
T ~ 10 10 yr). As it is apparent from eq. (|7^), t reac t is shorter for softer binaries. 
Hence, in order to contemplate hard binaries with an arbitrary orbital radius 
we have to set a x ~ 1.4 x 10 12 (Mq/M dc ) 2/3 km in eq. @. Correspondingly, 
the binary survival condition reads p*| 

'Afnr7\ 2 /M pc- 3 ~ ° 



1 (P) 



/pb < 8.2x10-°/cl S 
1 + 35/clS" 1 



'DC 

Mm 



P*(0) 



(0) \ / m y 5/3 



1 



(80) 



^Mopc" 3 , 

Although the presence of various dark cluster parameters prevents a clear - cut 
conclusion to be drawn from eq. (|80|), in the illustrative case Mdc — 1O 5 M 



and fcL ^ 0.5 eq. Q entails e.g. / PS <0.3 for ^(0) ~ 3 x 10 3 M Q pc" 3 . 
Thus, we infer that for realistic values of the parameters in question a sizable 
fraction of primordial binary brown dwarfs survives binary-binary encounters in 
the dark cluster cores F]. 

5 Thermal balance in halo molecular clouds 



As far as the energetics of halo molecular clouds is concerned, we can identify 
two main heat sources. One is energy deposition from background photons, 
while the other is of gravitational origin (coming from frictional hardening of 
primordial binary brown dwarfs) . Although we shall discuss them separately, it 
should be kept in mind that they both act at the same time. 

19 Application of the same argument to tidally - captured binaries shows that no depiction 
occurs in this way. 

20 If the initial value of fpB fails to satisfy condition (|sc|), primordial binaries start to be 
destroyed in binary-binary encounters until their fractional abundance gets reduced down to 
a value consistent with condition (po|). 
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5.1 Energy from background photons 

We proceed to estimate T m by momentarily neglecting gravitational effects. To 
this end, we need to know the heating rate (due to external sources) and the 
cooling rate (due to the molecules). In the galactic halo, the dominant heat 
source for molecular clouds is expected to be ionization from photons of the 
X-ray background, whose spectrum in the relevant range 1 keV < E < 25 keV 
(see below) can be parametrized in terms of the energy E (expressed in keV) 
as 1(E) = 8.5 E-° A cmT 2 s" 1 sr" 1 (O'Dea et al. ]1994| ). The ionization rate 
per Hi molecule (taking secondary ionization into account) is 

6c = 26 / 4irI(E)a x {E)dE s" 1 , (81) 

J E min 

where ax(E) — 2.6 x 10~ 22 E~ 8 / 3 cm 2 is the absorption cross-section in the 
above energy range for incoming X-rays on gas with interstellar composition 



(Morrison and McCammon [ 1983 1 ) . Actually, <jx(E) decreases faster as E in- 
creases when the gas metallicity is lower (see Fig.l in Morrison and McCammon 



[1983 1 ) . So, the above quoted expression for ax(E) is expected to yield un upper 
bound on the cross-section for X-rays on halo molecular clouds. The integration 
limits E min (below which X-rays are absorbed) and E max (above which X-rays 
go through the whole cloud without being absorbed) depend on the cloud col- 
umn density. Thus, taking as an orientation n m in the range 10 — 10 s cm -3 , we 
get 1.25 keV < E rnin < 7 keV and 10 keV < E max < 20 keV. It turns out that 
£x is rather insensitive to the upper limit, and we obtain £x = 2.2 x 10 -19 s _1 
for E m in = 1-25 keV and £ x = 5.6 x 10~ 21 s" 1 for E mm = 7 keV 0. Since 

21 It is straightforward to verify that the ionization rate £ cr due to the halo cosmic rays is 
less important in this context. In fact, assuming the ionization rate per Hi m olecule typical 
for disk vcosmic rays £o = 5 X 10 — 17 s — 1 (see e.g. van Dishoeck and Black ]19Sfi[) an d rescaling 
for the diffusion of the cosmic rays in the galactic halo (De Paolis et al! [ 1995b |), we infer 
£ cr - lO-M s-1. 
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each ionization process releases an energy of ~ 8 eV, the heating rate per Hi 
molecule T turns out to be 

3.5 x 1CT 32 erg s" 1 / 'H 2 < T < 1.5 X 1CT 30 erg s" 1 /-^ • (82) 

Unfortunately, the cooling rate for the halo clouds in question is not well 
known, owing to the lack of detailed information on their chemical composition. 
Nevertheless, by merely considering the cooling rate due to Hi as given by 



Goldsmith and Langer 1978], the equality between cooling and heating rate 
per molecule leads to T m ~ 10 K. More complete models for the cooling rate, 
which include the contribution from HD and heavy molecules, imply that the 
cooling efficiency is substantially enhanced and thus make it very plausible that 
halo molecular clouds clumped into dark clusters should have a temperature 



close to that of the CBR, namely T m ~ 3 K (see also Gerhard and Silk [ 1996 1 
for similar conclusions). 

5.2 Energy from primordial binary brown dwarfs 

As the analysis in Sect 4.3 shows, dynamical friction transfers a huge amount 
of energy from primordial binary brown dwarfs to molecular clouds, and so it 
looks compelling to investigate (at least) the gross features of the ensuing energy 
balance. 

Let us start by evaluating the energy acquired by molecular clouds in the 
process of frictional hardening. Recalling that the traversal time for a single 
cloud is given by eq. (|6l|), eq. (|58| ) entails that - after a binary with initial 
orbital radius a\ has crossed N clouds - its orbital radius gets reduced down to 



a-N+i 



r,n J V M DC J \ Ol / 



km . (83) 
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tfo S 9xlO»H(W^) (-) • (") 



Accordingly, we see that the orbital radius remains almost constant until N 
reaches the critical value 

V3 - X 3/2 

pc J V M Q / \ 01 

whereas it decreases afterwards. Because the energy acquired by the clouds is 
just the binding energy given up by primordial binaries, this information can be 
directly used to compute the energy AE C (N) gained by the N-th. cloud traversed 
by a binary whose initial orbital radius was eti. Manifestly, we have 

AE C (N) = ]-Gm 2 ( — -) . (85) 

I \a,N+l UN/ 

Thanks to eq. (|8Sj), a straightforward calculation shows that AE C (N) stays 
practically constant 

as long as N <Nq, while it subsequently decreases. So, the amount of energy 
transferred to a cloud is maximal during the early stages of hardening. Now, 
since the binding energy of a cloud is 

E c ~ 7.7 x 10 42 erg, (87) 

it can well happen that AE C > E c (depending on r m , M^c and a±), which 
means that the cloud would evaporate unless it manages to efficiently dispose 
of the excess energy. 

A deeper insight into this issue can be gained as follows (we focus the at- 
tention on the early stages of hardening, when the effect under consideration 
is more dramatic). Imagine that a spherical cloud is crossed by a primordial 
binary which moves along a straight line, and consider the cylinder A traced by 
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the binary inside the cloud (its volume being approximately 7ra 2 r m ). Hence, by 
n m ~ 62.2 (pc/r m ) 2 cm -2 (which follows from eq. (12)) the average number of 
molecules inside A turns out to be 

****** m 

Physically, the energy AE C gets first deposited within A in the form of heat. 
Neglecting thermal conductivity (more about this, later), the temperature inside 
A accordingly becomes 

ks being the Boltzmann constant. On account of eqs. © and ©, eq. (p9l) 
yields 

(90) 

which - in the illustrative case fcL — 0.5, p*(0) ~ 3 x 10 3 M S pc~ 3 and Mdc — 
1O 5 M -entails in turn 5.3 xl0 3 K <T A < 1.3xl0 4 K. As a consequence of the 
increased temperature, the molecules within A will radiate, thereby reducing 
the excess energy in the cloud. In order to see whether this mechanism actually 
prevents the cloud from evaporating, we notice that the characteristic time 
needed to accumulate the energy AE C inside A is just the traversal time t m . 
Therefore, this energy will be totally radiated away provided the cooling rate 
per molecule A exceeds the critical value Ao given by the equilibrium condition 

N A A t m ~ AE C . (91) 

Specifically, eq. Q) yields 

Ao ~ 10- 12 S- 1 (^j ^ erg s" 1 mob 1 , (92) 
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on account of eqs. (^lj), (^6|) and (|8|). Moreover, in the present case in which 
most of the molecules are H% the explicit form of A is (see e.g. O'Dea et al. 
1994, Neufeld et al. 1995) 

A ~ 3.8 x 1CT 31 erg s" 1 H^ 1 , (93) 

which - thanks to eq. (]89|) - becomes 



\ 4.35 

A~3.4x 10 21 S- 2 - 9 ( ™J ergs" 1 ^ 1 (94) 

(notice that A is almost independent of Mnc)- Now, from eqs. (|92| ) and (H) it 
follows that the condition A > Ao implies 

/ \ 0.35 

01 < 6x 10 11 E-°- 7 r^) km. (95) 

For a wide range of dark cluster and molecular cloud parameters it turns 
out that eq. ( psf ) is fulfilled for hard primordial binaries. 

Thus, we conclude that the energy given up by primordial binary brown 
dwarfs and temporarily acquired by molecular clouds is efficiently radiated away, 
so that the clouds do not get dissolved by frictional hardening. 

As a final comment, we stress that our estimate for Ta should be understood 
as an upper bound, since thermal conductivity has been neglected. In addition, 
the above analysis implicitly relies upon the assumption Ta < 10 4 K, which 
ensures the survival of H2. Actually, in spite of the fact that condition ( po| ) 
entails that this may well not be the case, our conclusion remains nevertheless 
true. For, higher temperatures would lead to the depletion of H2, which corre- 
spondingly gets replaced by atomic and possibly ionized hydrogen. As is well 
known, in either case the resulting cooling rate would exceed the one for H%, 
and so cooling would be even more efficient than estimated above (Bohringer & 
Hensler 1989). 
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6 Lyman-a absorption systems 

It is well known that Quasar Ly-a absorption lines provide a detailed informa- 
tion on the evolution of the gaseous component of galaxies (see e.g. Fukugita, 



Hogan and Peebles [1996] and references therein). These lines are seen for a neu- 
tral hydrogen column density Nhi ranging from ~ 3x 10 12 cm -2 (the detection 
threshold) to ~ 10 22 cm~ 2 . 

At the upper limit of this range (Nhi > 2 x 10 20 cm -2 ) the lines are classified 
as damped Ly-a systems (mostly associated with metal-rich objects) and it is 
generally believed that they are the thick progenitors of galactic disks. The HI 
distribution in damped Ly-a systems is usually flatter than the corresponding 
surface brightness of the optical disks and extends to much larger radii (up 
to ~ 40 kpc in giant galaxies). In the outer galactic regions Ly-a systems 
show sharp HI edges at a level of Nhi ~ 10 17 cm -2 , the so-called Lyman 
limit. Damped Ly-a systems are observed up to redshift z ~ 3.5 and the survey 
results suggest that the average mass of neutral hydrogen per absorption system 
decreases with time, in agreement with the hypothesis of gas consumption into 
stars (Lanzetta, Wolfe and Turnshek | 1995[ ). 



Within our picture, it is tempting to identify damped Ly-a systems with the 
PGC clouds in the inner galactic halo, where they undergo disruption. 

Below the Lyman limit (i.e. for 3 x 10 12 cm -2 < Nhi < 3 x 10 15 cm~ 2 , see 



e.g. Fig. 2 in Cristiani [ 1996 1 ) , the absorption lines are classified as Ly-a forest 
and are generally ascribed to a large number of intervening clouds (in some 
cases extending up to ~ 200 kpc from a central galaxy) along the line of sight to 
distant Quasars. Studies of Ly-a forest lines have made rapid progress recently 
and some observational trends are today firmly established. In particular: (i) 
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The evolution of the co-moving number density of systems per unit interval in 
redshift shows that the number of Ly-a forest clouds rapidly decreases with time 
for z > 2, while it is approximately constant for z < 2. (ii) The correlation 
between the thermal Doppler parameter b and the number of clouds can be fitted 
with a gaussian distribution of median < b > ~ 30 km s _1 and dispersion ~ 10 
km s -1 , corresponding to a temperature of a few 10 4 K. 

Within our model, it looks natural to identify Ly-a forest clouds with the 
molecular clouds clumped into dark clusters located in the outer halo. Indeed, 
we expect that the clouds contain in their external layers an increasing fraction 
of HI gas and that the outer regions are even ionized, due to the incoming UV 
radiation Observe that so far we have been assuming that halo molecular 
clouds have a constant temperature. This was a good approximation as far 
as the previous analysis was concerned, but would be too poor in the present 
discussion. A HI column density of ~ 10 14 cm -2 - corresponding to a layer 
of about 10~ 6 pc - is sufficient to shield the incoming radiation. Remarkably 
enough, Nhi ~ 10 14 cm~ 2 corresponds to the average value of the observed 



Lyman-a forest distribution (Cristiani [1996 1). Since we expect the UV flux to 
decrease as the galactocentric distance increases, clouds which lie at a larger 
distance may thus have a smaller HI column density, whereas clouds closer 
to the galactic disk may have a higher column density. This fact explains the 
observed distribution for the column density. Moreover, the number of clouds is 
expected to decrease with time initially, due to both MACHO and H2 formation, 
in this way explaining the observed evolution of Ly-a forest clouds according 

22 Observe that so far we have been assuming that halo molecular clouds have a constant 
temperature. This was a good approximation as far as the previous analysis was concerned, 
but would be too poor in the present discussion 
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to the above-mentioned point (i). As a final comment, we mention that very 
recently Rottgering, Miley and Van Ojik 199(| , by considering the filling factors 



and the physical parameters derived from their Ly-a forest observations, pointed 
out that galactic halos may contain ~ 10 9 Mq of neutral hydrogen gas and are 
typically composed of ~ 10 12 clouds, each of size ~ 40 light-days. It looks 
remarkable that these are the typical parameters of molecular clouds dealt with 
in this paper. 

A thorough quantitative analysis requires, however, further investigations, 
which are beyond the scope of the present paper. 

7 Conclusions 

Looking back at what we have done, it seems to us fair to say that present- 
day results of microlensing experiments towards the LMC are amenable of a 
very simple explanation in terms of what was repeatedly proposed as a natural 
candidate for baryonic dark matter, namely brown dwarfs. Once this idea is 
accepted, a few almost obvious steps just follow. First, given the fact that 
ordinary halo stars form in (globular) clusters, it seems more likely that brown 
dwarfs as well form in clusters rather than in isolation (this circumstance has 
been repeatedly recognized in the last few years). Of course, whether brown 
dwarfs are still clumped into dark clusters today is a different and nontrivial 
question - we have seen that core collapse can liberate a considerable fraction of 
brown dwarfs from the less massive dark clusters. Second, much in the same way 
as it happens for ordinary stars, a consistent fraction of binary brown dwarfs 
(up to 50% in mass) is expected to form. Third, since no known star-formation 
mechanism is very efficient, it is natural to imagine that a substantial amount 
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of primordial gas is left over. Although we cannot be sure about the subsequent 
fate of this gas - which should be mostly cold if 2 - it is likely to remain confined 
within the dark clusters, since stellar winds are absent. So, also cold H2 self- 
gravitating clouds should presumably be clumped into dark clusters, along with 
some residual diffuse gas. 

We have shown that - within the considered model - not too hard primordial 
binary brown dwarfs turn out to have an orbital radius which is typically of the 
order of the Einstein radius for microlensing towards the LMC. Therefore, they 
are not so easily resolvable in the microlensing experiments performed so far, but 
we argue that they can be resolved in future microlensing experiments with a 
more accurate photometric observation ^|, the signature being small deviations 
from standard microlensing light curves (Dominik 1996fl ). Notice that such a 
procedure complements the detection strategy for binaries suggested by Gaudi 



and Gould |1997 |. 

Many physical processes certainly occur in the dark clusters, and the (likely) 
presence of a large amount of gas makes them even more difficult to understand 
than globular clusters. Moreover, the lack of any observational information 
makes any attempt to figure out the physics of dark clusters almost impossible, 
were it not for some remarkable structural analogies with globular clusters that 
the model - upon which our discussion is based - naturally suggests. What 
should anyway be clear is that a host of different phenomena can happen. Some 
of them have intentionally been neglected here, in order not to make the present 
paper either excessively long or too speculative. Yet, many others are likely to 

23 See e.g. the ongoing experiments by the GMAN and PLANET collaborations (Proceed- 
ings of the Second International Workshop on Gravitational Microlensing Surveys, Orsay, 
1996). 
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occur, that we have not even been able to imagine. 
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